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INTRODUCTION 


The  purpose  of  this  paper  is  to  provide  a  qualitative  assessment  of 
total  electron  content  (TEC)  forecasting  at  the  Air  Force  Global  Weather 
Central  (AFGWC).  It  is  based  on  five  years  of  forecasting  experience 
encompassing  the  extremes  of  TEC  variation  between  solar  minimum  and  solar 
maximum.  TEC  products  are  generally  of  two  types:  an  area  forecast  or  a 
point  forecast.  Area  forecasts  are  for  geographical  areas  on  the  order  of  ten 
degrees  of  latitude  and  longitude  or  larger.  In  general,  the  desired  goal  is 
to  predict  or  specify  TEC  withinj^  5  X  10^  electrons/m^.  Throughout  this 
paper  the  various  aspects  of  the  prediction  and  specification  process  will  be 
examined.  A  detailed  discussion  of  TEC  modeling  techniques  is  beyond  the 
scope  of  this  paper.  Instead  the  difficulties  and  subjectivity  of  the 
forecast  process  will  be  emphasized. 

OBSERVATIONAL  DATA 


The  available  data  set  includes  hourly  TEC  observations  from  12 
polarimeter  sites  that  monitor  VHF  beacons  on  geosychronous  satellites.  Seven 
of  these  locations  are  24  hour  per  day  reporting  sites.  The  remaining  sites 
operate  during  normal  duty  hours  and  report  the  remaining  observations  in  a 
daily  summary.  This  "network"  of  TEC  observing  sites  is  complemented  by 
ground  based  vertical  incidence  (V.I.)  sounders.  Six  V.I.  sounder  sites 
transmit  data  on  a  24  hour  per  day  basis.  Summary  data  are  received  between 
12-72  hours  after  the  observation  times  from  an  additional  25-30  stations. 

TEC  observations  have  proven  to  be  the  key  input  parameter  to  the 
forecast.  The  TEC  observations  have  a  number  of  problems,  foremost  of  which 
is  the  lack  of  sufficient  observing  sites.  For  area  forecasts  one  would 
ideally  like  to  have  the  forecast  area  boxed  in  by  observing  sites.  For  areas 
that  cover  a  large  latitudinal  variation,  a  chain  of  stations  along  the  same 
meridian  is  desirable.  A  good  example  of  this  is  the  chain  of  polarimeter 
sites  along  the  North  American  east  coast  consisting  of:  Goose  Bay,  Labrador; 
Sagamore  Hill,  Massachusetts;  Patrick  AFB,  Florida;  and  Ramey  AFB,  Puerto 
Rico.  This  chain  can  effectively  monitor  the  high  and  middle  latitude 
ionosphere  and  lacks  only  an  equatorial  station.  Such  chains  of  observing 
sites  near  the  same  meridian  record  the  ionosphere  near  the  same  local  time 
and  greatly  aid  the  forecaster  in  identifying  synoptic  patterns  that  exist  in 
a  local  time  frame.  Three  or  four  such  meridional  chains  are  needed  to 
monitor  TEC  on  a  worldwide  basis.  For  forecasts  in  areas  of  high  ionospheric 
variability  such  as  the  auroral  zone  or  equatorial  region,  additional  TEC 
sites  are  needed  as  well  as  more  frequent  observations.  Observations  from 
these  areas  should  routinely  be  transmitted  every  fifteen  minutes.  During 
disturbed  conditions,  fifteen  minute  observations  should  also  be  received  from 
middle  latitude  stations. 

Data  receipt  and  reliability  are  other  problems  the  forecaster  has  to 
contend  with.  Overseas  locations  often  experience  communication  outages  or 
delays.  This,  combined  with  the  small  number  of  observing  sites,  can  lead  to 
a  very  sparse  data  base  at  times.  The  quality  of  the  observations  are 
occasionally  suspect  or  in  error,  because  the  observer  personnel  lack  both  the 
knowledge  of  ionospheric  variabi 1 ities  and  training  needed  to  interpret  the 
data.  This  is  particularly  true  during  disturbed  ionospheric  conditions  when 
rapid  variations  occur  that  make  TEC  measurements  difficult  for  the  best  of 
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observers.  The  result  of  these  erroneous  observations  is  that  TEC 
observations  from  a  given  station  may  be  in  error  by  5-10  X  10^ 
electrons/m^  for  periods  ranging  from  several  hours  to  several  days.  Such 
values  exceed  the  desired  forecast  goals,  and  the  erroneous  data  is  difficult 
to  identify  due  to  the  lack  of  corroborating  ionospheric  data. 

Several  of  the  polarimeter  sites  are  completely  automated.  While 
this  has  the  advantage  of  24  hour  per  day  observations,  the  software 
occasionally  has  difficulties  in  keeping  track  of  the  TEC  values  during 
disturbances  and  the  sunrise  transition  period.  With  unmanned  sites,  the 
correction  of  such  problems  is  even  more  difficult.  An  additional  problem, 
which  is  becoming  increasingly  more  critical,  is  the  declining  availability  of 
VHF  beacons  on  geosychronous  satellites.  Not  only  does  this  limit  the 
geographic  extent  of  the  observing  network,  but  the  operation  and  maintenance 
of  the  satellites  and  their  beacons  are  not  under  the  control  of  the  TEC 
monitoring  network.  Add  to  this  the  lack  of  standardized  observing  equipment 
and  centralized  control,  and  the  problems  associated  with  the  operation  and 
maintenance  become  more  complex.  From  the  forecaster's  standpoint,  this 
results  in  an  excessive  amount  of  time  spent  in  data  acquisition  and 
validation. 

FORECAST  MODELS 


There  are  several  approaches  to  input  data  into  the  AFGWC  4-D 
ionospheric  model.  TEC  data  can  be  used  as  the  primary  input  parameters,  or 
it  can  be  combined  with  foF2  data  to  produce  TEC  specification  fields. 

Another  option  is  to  use  only  foF2  data.  A  fourth  option  is  to  use  the 
Institute  of  Telecommunications  Service  (ITS)  data  fields.  It  should  be  noted 
that  all  of  the  ionospheric  models  in  use  are  specification  models  that  have 
no  physical  dynamics  built  into  them.  Background  information  concerning 
ionospheric  modeling  efforts  at  AFGWC  is  contained  in  a  paper  by  Tascione  et 
al,  1979  (see  references).  The  TEC  forecast  products  are  a  result  of 
forecaster  modification  of  the  model  or  model  output  based  on  experience. 

On  a  given  day  any  one  or  more  of  the  various  input  parameters  may 
produce  the  best  TEC  specification  field.  Over  the  long  term;  however,  the 
option  using  foF2  and  TEC  data  produces  the  best  TEC  specification.  This 
method  uses  considerable  computer  time  and  produces  only  a  slightly  better  TEC 
field  than  the  use  of  TEC  alone  as  an  input  parameter.  The  use  of  foF2 
observations  produces  good  TEC  specification  on  some  days  but  lacks  the 
consistency  for  use  on  a  day-to-day  basis.  The  ITS  fields  also  suffer  from  a 
similar  lack  of  consistency. 

TEC  fields  produced  from  the  ITS  model  are  not  viable  for  routine 
daily  forecasts.  Predictions/specifications  from  the  ITS  model  are  often 
mistaken  for  climatology.  There  is  no  TEC  climatology  except  for  the  data 
currently  being  collected.  The  basic  input  for  the  ITS  model  is  a  sunspot 
number  (SSN).  The  best  available  sunspot  number  for  ionospheric  purposes  is 
the  effective  sunspot  number  (SSNeff)  produced  from  actual  foF2 
observations.  Table  l  shows  the  extreme  variation  in  daily  effective  sunspot 
numbers  calculated  through  the  month.  The  variations  between  low  and  high 
daily  values  within  a  given  month  range  from  50-60  to  over  100.  These 
effective  sunspot  numbers  are  the  best  fit  of  an  ITS  field  to  the  observed 
foF2  values  for  a  given  day.  The  use  of  an  ITS  field  specified  on  other  than 
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a  daily  basis  could  result  in  extremely  poor  TEC  fields  depending  upon  the 
particular  ionospheric  conditions  on  a  given  day. 


TABLE  I 


Daily  Effective  Sunspot  Number 


MONTH 

LOW 

HIGH 

Jul  1978 

72 

142 

Aug 

60 

119 

Sep 

54 

117 

Oc  t 

84 

137 

Nov 

84 

132 

Dec  1978 

89 

118 

Jan  1979 

112 

140 

Feb 

142 

193 

Mar 

114 

200 

Apr 

88 

200 

May 

99 

181 

Jun 

107 

172 

Jul 

112 

145 

Sep 

121 

140 

Oct 

122 

184 

Nov 

131 

181 

Dec  1979 

121 

177 

Jan  1980 

124 

158 

Feb 

123 

160 

To  illustrate  some  of  the  TEC  variations  consider  the  following.  For 
February  1980,  a  TEC  field  generated  from  ITS  fields  with  sunspot  numbers  of 
120,  140,  and  175  at  a  time  of  1100Z  would  produce  the  following  TEC 
variations . 


TABLE  2 


ITS  TEC  (10^  electrons/m^ ) 


LOCATION 

SSN; 

120 

140 

175 

Athens 

48 

57 

76 

Osan 

26 

32 

44 

Taiwan 

68 

81 

105 

Observed  values  at 

11002  during 

Feb  1980 

we  re : 

LOCATION 

Minimum 

Maximum 

Athens 

35 

72 

Osan 

14 

34 

Taiwan 

51 

125 
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The  observed  values  above  are  not  the  extremes  associated  with  large 
disturbances  but  typical  of  daily  variations  within  a  given  month.  If,  for 
example,  TEC  predictions  were  issued  on  ;t  monthly  basis  usinft  the  ITS  fields, 
large  differences  between  predicted  and  observed  TEC  values  could  be  expected 
during  undisturbed  conditions.  Use  of  monthly  predictions  also  neglects  the 
following : 


1.  The  day-to-day  variablity  of  the  equatorial  and  auroral  zone 
regions  is  not  considered. 

2.  Quiet  day  TEC  values  and  patterns  are  often  quite  different 
between  the  Eastern  and  Western  Hemispheres. 

3.  Seasonal  changes  in  the  ionospheric  characteristics  can 
occur  at  different  times  then  the  calendar  season  change. 

During  disturbed  ionospheric  conditions,  ITS  fields  are  totally 
unrepresentative.  TEC  variations  can  be  extreme.  For  example,  during  a 
storm,  values  at  Osan,  Korea  can  be  30  TEC  units  (1  TEC  unit  =  1.0  X  10 ^ 
electrons/m^)  below  the  previous  quiet  day  values.  At  the  same  time,  TEC 
values  at  Athens,  Greece  can  be  30  TEC  units  above  the  previous  quiet  day's 
value.  On  one  occasion,  nighttime  values  at  Palehua,  Hawaii  jumped  nearly  25 
TEC  units  (typical  values  are  5-8  TEC  units)  for  several  hours.  TEC  values 
at  Taiwan  have  changed  from  values  of  >  150  TEC  units  on  a  quiet  day  to  values 
of  <  50  on  the  following  disturbed  day. 

QUIET  DAY  TEC  FORECASTS 

The  ionosphere  can  show  remarkable  variations  even  during  undisturbed 
conditions.  The  quiet  day  variations  sometimes  match  storm  variations.  That 
is,  sometimes  large  increases  or  decreases  are  noted  for  a  day  or  more  with  no 
significant  geomagnetic  disturbances  in  progress.  Generally,  the  rapid 
variations  often  observed  during  geomagnetic  storms  are  lacking  but  the 
magnitude  of  the  TEC  change  can  be  comparable.  Large  gradients  of  TEC  are 
also  apparent  in  certain  areas.  For  example,  on  some  days  Palehua,  Hawaii  TEC 
values  are  similar  to  middle  latitude  stations.  On  other  days  it  resembles  an 
equatorial  ionospheric  station  such  as  Taiwan.  It  should  also  be  noted  that 
despite  the  relatively  small  geographical  separation,  daytime  TEC  values  at 
Taiwan  are  nearly  double  those  of  Osan,  Korea.  Nighttime  values  show  even 
greater  differences. 

The  approach  to  TEC  forecasting  during  quiet  conditions  is  to  have 
the  model  build  a  specification  field  based  on  persistence  or  five  day  mean 
TEC  values.  The  forecaster  lias  some  flexibility  in  use  of  the  model  so  that 
unrepresentative  TEC  for  various  stations  or  days  can  be  excluded  from  the 
desired  specification  field.  By  monitoring  and  displaying  the  TEC  data  on  a 
routine  basis  the  forecaster  attempts  to  identify  synoptic  patterns  and  trends 
for  application  to  the  specification  fields  produced  by  the  model.  The 
forecaster  can  then  apply  the  changes  necessary  to  modify  the  model  output  or 
the  TEC  forecast  itself.  In  the  case  of  TEC  forecasts  over  large  areas  and/or 
time  periods,  the  forecaster  can  make  changes  by  area  or  local  time 
application  of  weighting  factors.  These  factors  then  produce  the  desired 
increases  or  decreases  in  predicted  TEC  values. 
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DISTURBED  DAY  TEC  FORECASTS 


During  disturbed  ionospheric  conditions  the  forecast  problem  is  much 
more  difficult.  Not  all  geomagnetic  disturbances  produce  observed  large  scale 
TEC  variations.  At  other  times  similar  geomagnetic  disturbances  produce 
completely  different  TEC  responses-  There  appears  to  be  a  marginal 
correlation  between  the  magnitude  of  the  storm  and  the  magnitude  of  the  TEC 
variations.  The  sparsity  of  data  precludes  early  detection  and  continual 
monitoring  of  the  disturbance.  Figures  1  through  6  illustrate  the  type  of  TFC 
variation  that  is  possible  during  a  disturbance.  The  forecast  process  is 
hampered  by  the  lack  of  a  dynamic  prediction  model  and  essentially  no 
prediction  tools  for  the  forecaster.  Storm  morphology  in  the  research 
literature,  until  recently,  has  concentrated  on  the  middle  latitude  and  is 
generally  based  on  a  single  station.  As  a  general  rule  the  forecaster  will 
look  for  TEC  increases  (positive  phase)  of  10  to  10  percent  above  the  quiet 
day  values  early  in  the  disturbance  followed  by  a  sharp  decrease  (negative 
phase)  to  below  quiet  day  values.  The  enhancement  figure  of  30  percent  above 
the  quiet  day  value  appears  o  be  a  common  value  for  middle  latitudes.  The 
timing  of  positive  and  negative  phases  is  very  difficult.  On  a  global  basis, 
one  hemisphere  may  be  in  one  phase  while  the  other  hemisphere  experiences  the 
opposite  phase.  Often  the  forecast  is  out  of  phase  with  the  observed  values. 

FORECAST  VERIFICATION 


The  verification  of  point  forecasts  in  areas  where  observational  data 
exists  is  straightforward  ana  can  be  done  in  a  variety  of  ways.  In  areas 
where  there  are  no  TEC  observations,  feedback  from  the  user  of  the  TEC 
forecast  is  the  only  mechanism  for  evaluating  the  forecast.  Since  roost  users 
are  interested  in  TEC  forecasts  over  an  area  rather  than  a  single  point,  the 
use  of  observations  from  a  single  station  is  very  risky  if  the  area  extends 
very  far  in  either  latitude  or  longitude.  Some  attempts  have  been  made  to. use 
observations  from  a  few  stations  as  verification  data  representative  of  a 
given  latitude  regime  on  a  worldwide  basis.  For  example,  a  TEC  forecast  for 
the  middle  of  the  Pacific  Ocean  would  be  verified  for  the  same  local  time 
using  a  middle  latitude  station  on  the  east  coast  of  the  IT. S.  These  schemes 
neglect  the  fact  that  the  forecaster  modifications  are  for  the  forecast  area 
not  Che  area  of  the  verification  data.  Such  efforts  are  analogous  to  using 
Omaha,  Nebraska  weather  observations  to  verify  a  weather  forecast  for  Los 
Angeles,  California.  Until  a  hi gh  density  network  of  observations  is 
available,  the  only  practical  verification  scheme  is  some  measure  of  the 
forecast's  value  by  the  user. 

SOLAR  CYCLE  EFFECTS  ON  TEC  FORECASTING 


A  few  comments  concerning  the  variation  ot  TEC  during  the  solar  cycle 
and  its  effect  on  TEC  forecasting  should  be  made.  Figure  7  illustrates  the 
typical  differences  of  TEC  values  during  solar  minimum  and  maximum 
conditions.  Maximum  daily  TFC  values  at  the  middle  latitude  during  solar 
minimum  are  on  the  order  of  15-2  3  X  101*5  electrons/-.!^  compared  to  values 
of  60-70  X  lO1^  electrons/m--  d  uring  solar  maximum.  Even  more  extreme 
variations  are  exhibited  ,n  the  equatorial  region  where  peak  daily  values  of 
50  X  10^  e  lectrons/ni?  during  solar  minimum  changed  to  values  of  150  X 
10*6  electrons/m^  dining  solar  maximum.  Nighttime  values  of  TEC.  remain 
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Notable  exceptions  to  those  values  may  occur  during  disturbed 
ionospheric  conditions  or  other  unusual  solar  conditions,  but  the  nighttime 
values  do  not  approach  the  daytime  values. 

The  large  daily  variations  in  TEC  during  solar  maximum  also 
complicate  the  TEC  observing  process.  The  large  number  of  phase  changes  that 
have  to  be  determined  using  the  Farady  rotation  technique  adds  to  the 
likelihood  of  erroneous  observations.  This,  in  turn,  adds  an  additional 
workload  to  the  forecaster's  effort  to  have  an  error-free  database. 

SUMMARY 

TEC  forecasting  at  AFGWC  has  been  limited  by  a  number  of 
difficulties.  One  of  the  primary  problems  is  observational  data.  There  are 
insufficient  data  for  model  input  or  for  forecast  verification.  The  existing 
polarimeter  network  is  hampered  by  training  problems  and  lack  of  centralized 
control.  This  produces  an  additional  forecaster  workload  in  maintaining  a  TEC 
data  base.  To  eliminate  this,  a  future  observing  network  should  consist  of 
fully  automated,  dual  frequency  TEC  measurements. 

Because  of  data  problems,  the  capabilities  of  AFGWC  are  essentially 
untested  for  area  forecasts.  Significant  advances  in  the  specification  model 
have  been  made,  but  its  success  in  an  operational  environment  is  unknown 
because  of  insufficient  observations  and/or  customer  feedback.  The  solar 
cycle  variations  in  TEC  make  it  doubtful  that  the  desired  accuracies  could  be 
achieved  during  solar  maximum.  With  sufficient  observations,  the  desired 
prediction  accuracies  appear  achievable  during  solar  minimum.  Despite  this, 
daily  TEC  predictions  using  the  existing  data  set  are  a  significant 
improvement  over  models  that  are  updated  only  on  a  monthly  basis. 

For  point  forecasts,  prediction  goals  can  be  met  during  undisturbed 
conditions  if  representative  TEC  observations  are  available  near  the  area  of 
interest.  However,  day  to  day  variations  during  solar  maximum  will  test  the 
forecaster's  skill.  Forecaster  skill  can  only  be  obtained  through  a  thorough 
knowledge  of  ionospheric  morphology  and  ionospheric  forecasting  experience. 
Because  the  physical  processes  are  more  complex  than  conventional  meteorology 
and  predictive  techniques  virtually  non-existent,  forecaster  experience  in 
identifying  synoptic  patterns  and  trends  in  the  ionosphere  is  critical  to  the 
success  of  the  forecast.  TEC  forecasting  during  disturbed  conditions  shows 
little  or  no  skill.  Observational  data  problems  and  forecaster  inexperience 
contribute  to  this,  but  inadequate  resources  have  been  devoted  to  this 
difficult  aspect  of  the  forecast  problem.  Some  common  features  appear  at 
various  locations,  but  their  relationship  to  other  geographical  parameters  is 
not  clear. 
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